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APPROXIMATOR APPARATUS, AND ASSOCIATED METHOD, FOR 
APPROXIMATING AN N-DIMENSIONAL EUCLIDEAN NORM 

The present invention relates generally to a manner by which to approximate an 
5 N-dimensional Euclidean norm. More particularly, the present invention relates to 
approximator apparatus, and an associated method, by which to approximate the value 
of an N-dimensional Euclidean norm of a data symbol defined in terms of N dimensions. 

When implemented to approximate the Euclidean norm of a two dimensional data 
symbol, increased accuracy of the approximation relative to the accuracy of 
10 approximation permitted by conventional approximation apparatus, is achievable. When 
implemented to approximate the Euclidean norm of a three dimensional data symbol, the 
approximator apparatus is of significantly reduced complexity relative to conventional 
approximation apparatus while maintaining high levels of accuracy. And, 
approximations of Euclidean norms of data symbols of higher level dimensions are also 
15 possible, made without necessitating significantly increased circuitry or processing 
complexity. 

The approximator apparatus is advantageously implemented in any 
communication device, such as a modem forming part of a sending or a receiving station 
of a CDMA (code-division, multiple-access) cellular communication system in which 
20 Euclidean norms are required to be formed to facilitate processing of data symbols of 
data sequences that are to be communicated during operation of the communication 
device pursuant to effectuation of a communication service, 



BACKGROUND OF THE INVENTION 

25 A communication system is formed, at a minimum, of a sending station and a 

receiving station, interconnected by a communication channel. Data that is to be 
communicated between the sending and the receiving stations is converted, by the 
sending station, into a form to permit its communication upon the communication 
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channel. And, the receiving station operates to recover the informational content of the 
data, representations of which are detected at the receiving station. 

Many varied types of communication systems have been developed and deployed 
to permit effectuation of many types of communication services. As technological 

5 v advancements permit, improvements to existing such communication systems are made, 
and new types of communication systems are developed and deployed. 

Amongst the technological advancements that have been implemented in 
communication systems have been advancements in digital communication techniques. 
Digital communication techniques, when implemented in a communication system, 

10 provide the capability of the communication system more efficiently to utilize the 

communication resources allocated thereto. By increasing the efficiency by which the 
communication resources of the communication system are utilized, the communication 
capacity of the communication system is increased. For instance, the bandwidths of the 
communication channel extending between the sending and receiving stations of a 

15 communication system are sometimes of limited allocations. That is to say, the 
communication capacity of such a communication system is constrained by the 
bandwidth available upon the communication channel upon which to communicate data. 

A radio communication system is an exemplary type of communication system. 
In a radio communication system, the communication channel that extends between the 

20 sending and receiving station is defined upon a radio link, i.e., a portion of the 

electromagnetic spectrum. Because data is communicated on radio channels, the need to 
utilize a wireline, conventionally required in a wireline communication system, upon 
which to define channels is obviated. Radio communication systems permit for the 
effectuation of communication services when formation of a wireline connection 

25 between the sending and receiving stations is inconvenient or impractical. Additionally, 
a radio communication system is implementable as a mobile communication system in 
which one, or both, of the sending and receiving stations is permitted mobility. 

A cellular communication system is a type of radio communication system that 
forms a mobile communication system. In a cellular communication system, mobile 

30 stations communicate by way of radio channels defined upon a radio air interface with 
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network based communication stations. Telephonic communication of data, both of 
voice and of non-voice data, is effectuable pursuant to conventional cellular 
communication systems. With the implementation of digital communication techniques 
in cellular communication systems, data-intensive communication services are 
5 effectuable. 

To facilitate communication of data in such digital cellular communication 
systems, as well as in other digital communications, data coding schemes have been 
devised. For instance, channel encoders are regularly utilized that assign symbols to 
binary representations of data that is to be communicated pursuant to the communication 

10 service. Corresponding channel decoders operate in reverse manners to convert detected 
symbol values into binary form. In some existing cellular communication systems, for 
instance, QPSK (quadrature phase shift keying) modulation is utilized in which symbols 
are assigned during channel encoding and modulation operations in which symbols are 
selected from an allowable, QPSK symbol set. Each symbol is defined in terms of an I- 

15 component and a Q-component. The QPSK scheme is a two-dimensional scheme. Other 
two-dimensional schemes are known and used in some other communication systems. 
By extension, three dimensional and higher dimensional schemes are also known and are 
amenable for implementation to facilitate effectuation of communication of data 
pursuant to a communication service. 

20 Communications devices, such as sending stations and receiving stations, of a 

digital communication system, are sometimes required to carry out computationally 
intensive operations. The need to carry out the computationally intensive operations 
limits the speed at which the communication operations can be carried out and requires 
complex, and therefore relatively costly circuitry, to carry out such operations. For 

25 instance, the operations that are sometimes required to be carried out upon the data 
include the calculation of Euclidean norms. Calculation of Euclidean norms requires 
both squaring and square root operations. Euclidean norm calculations are required to 
be formed, for instance, at modem devices that form parts of the communication stations 
of a cellular communication system. A modem forming part of a communication device, 

30 e.g., of a CDMA based (code-division, multiple-access-based) cellular communication 
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system, must perform Euclidean norm operations pursuant to operations of its various 
functions, such as searching, dll, transmitting, etc. functions. 

Calculating the Euclidean norms in a two dimensional scheme is difficult, 
particularly when such operations are required to be performed quickly and successively 
5 upon successive data symbols of data that is communicated to effectuate a 

communication service. Calculation of Euclidean norms in higher dimensional schemes 
become prohibitively complex, particularly at high data rates. 

Approximation techniques by which to approximate Euclidean norms are 
therefore sometimes used in existing systems and proposed for systems that shall utilize 

10 multidimensional techniques. Approximation techniques are also sometimes used in 
systems that utilize two-dimensional techniques. To be effective, any approximation 
technique must be accurate while also being of reduced computational complexity. 
While various approximation techniques are known to approximate two dimensional 
Euclidean norms, their application to higher dimensional Euclidean norm 

15 approximations also becomes computationally complex. As approximation techniques 
are presently not generally available for approximating quickly and accurately higher 
dimensional Euclidean norms, modem devices, as well as other devices, that are to be 
utilized to operate upon data that is represented in multiple dimensions either 
inaccurately provide approximations of Euclidean norms or are prohibitively complex 

20 and slow. 

Accordingly, an improved manner by which to approximate Euclidean norms that 
is both accurate and computationally noncomplex is required. 

It is in light of this background information related to digital communication 
systems, such as cellular communication systems, that the significant improvements of 
25 the present invention have evolved. 

SUMMARY OF THE INVENTION 

The present invention, accordingly, advantageously provides apparatus, and an 
associated method, by which to approximate an N-dimensional Euclidean norm. 
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Through operation of an embodiment of the present invention, a manner is 
provided by which to approximate the value of an N-dimensional Euclidean norm of a 
data symbol defined in terms of N dimensions. 

An approximation of the Euclidean norm is made in a manner that accurately 
5 represents the Euclidean norm while not requiring significant levels of computational 
complexity to form the approximations. 

When the Euclidean norm of a two dimensional data symbol is to be 
approximated, an approximation is formed that is of increased accuracy relative to the 
accuracy of an approximation made by a conventional approximation technique. 
10 When the Euclidean norm of a three dimensional symbol is to be formed, an 

approximation that is both accurate and of significantly reduced complexity relative to 
conventional approximation techniques is formed. 

And, when the Euclidean norm of a higher level, N-dimensional symbol is to be 
formed, an accurate approximation is obtainable without a corresponding significant 
15 increase in complexity of the computations required to form the approximation. When a 
N-dimensional Euclidean norm is to be approximated, N-l iterations of a two 
dimensional approximation are performed in which successive ones of the iterations 
approximate, iteratively, two dimensional Euclidean norms, successively utilizing 
successive ones of the component portions of the N-dimensional data symbol. 
20 Approximation of an N-dimensional Euclidean norm is thereby made merely through 
successive iterations of a two dimensional Euclidean norm approximation. Coordinates 
associated with a prior iteration of a two dimensional Euclidean norm are used for a 
following calculation of a two dimensional Euclidean norm approximation. 

In one aspect of the present invention, a first pair, or other set, of component 
25 portions of a data symbol is provided to a transformer. The pair, or other set, of 

component portions of the data symbol define an angular value. And, the transformer 
operates to transform, if appropriate, the values provided thereto into an angular value 
within a selected angular range, e.g., between 0 and pi/4 radians. 

And, further, the selected angular range is divided into a first angular subrange 
30 and a second angular subrange. The angular value into which the pair, or other set, of 
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component portions of the data symbol are transformed falls within one or the other of 
the angular subranges. The angular subrange into which the angular value defined by 
the transformed values falls is determinative of the manner by which the approximation 
of the Euclidean norm is made. If the angular value falls within the first angular 
5 subrange, the Euclidean norm is calculated in the first manner; if the angular value falls 
within the second angular subrange, the Euclidean norm is approximated in a second 
manner. 

In one implementation, the Euclidean norm approximations are performed at a 
modem device of a communication station operable in a CDMA cellular communication 

10 system. During operation of the modem device, successive data symbols forming data 
that is communicated to effectuate a communication service are applied to the modem 
device. Amongst the operations that are performed upon the data symbols are 
approximations of Euclidean norms. The approximator operates to form an 
approximation of the Euclidean norm of the data symbol, responsive to the values 

15 provided thereto. If a multi dimensional Euclidean norm is to be obtained, the 
approximation is formed pursuant to an iterative process in which successive two 
dimensional Euclidean norm approximation values are calculated such that after a 
desired iteration, the multi dimensional Euclidean norm approximation value is 
obtained. 

20 Because the approximation that is formed is both accurate and obtained without 

the need to perform complex computations, the approximator and its associated method 
of approximating is advantageously utilized when operations are to be performed upon 
large numbers of input data, such as data symbols that are to be communicated pursuant 
to a data communication service. When a two-dimensional Euclidean norm is to be 

25 approximated, improved accuracy relative to conventional approximation techniques is 
provided. And, higher level Euclidean norm approximations are formable while 
conventional techniques are computationally prohibitively high. 

In these and other aspects, therefore, an approximator, and an associated method, 
is provided for a communication device. The approximator forms an approximation of a 

30 Euclidean norm of a data symbol defined in terms of a first component portion and at 
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least a second component portion. A detector is adapted to receive indications of the 
first and at least second component portions of the data symbol. A transformer is 
coupled to the detector. The transformer transforms the component portion into a first 
transformed value and transduces the second component portion into a second 
5 transformed value. The first and second transformed values, respectively, formed by 
the transformer are selected such that a geometric argument defined therefrom is within 
a selected angular range. An estimator is adapted to receive the first and second 
transformed values. The estimator estimates the Euclidean norm of the data symbol. 
The Euclidean norm is estimated to be a summation of a first product formed from the 
10 first transformed value and a second product value formed from the second transformed 
value. 

A more complete appreciation of the present invention and the scope thereof can 
be obtained from the accompanying drawings that are briefly summarized below, the 
following detailed description of the presently-preferred embodiments of the invention, 
15 and the appended claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 illustrates a functional block diagram of a radio communication system 
in which an embodiment of the present invention is operable. 

Figure 2 illustrates a graphical representation of angular subregions defined 
20 pursuant to an embodiment of the present invention by which to approximate an N- 
dimensional Euclidean norm. 

Figure 3 illustrates graphical representations of error curves that indicate 
maximum error percentages of approximations generated during operation of an 
embodiment of the present invention. 
25 Figure 4 illustrates a graphical representation of a manner by which a Euclidean 

norm, of any of N dimensions is approximated pursuant to operation of an embodiment 
of the present invention. 

Figure 5 illustrates a method flow diagram of the method of operation of an 
embodiment of the present invention. 
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DETAILED DESCRIPTION 

Referring first to Figure 1, a radio communication system, shown generally at 10, 
provides for radio communications between mobile stations, of which the mobile station 
12 is representative, and a network part. In the exemplary implementation, the 
5 communication system forms a cellular communication system operable generally 
pursuant to an operating specification that utilizes code division, multiple-access 
(CDMA) communication techniques. The communication system is also representative 
of other types of cellular communication systems. And, more generally, the 
communication system 10 is representative of any of various other types of digital 

10 communication systems in which digital communication stations are interconnected by 
way of a communication channel. Accordingly, while the following description of 
exemplary operation of the present invention shall be described with respect to its 
implementation in which the communication system forms a CDMA-based, cellular 
communication system, the teachings of the present invention are analogously also 

15 applicable with respect to its implementation in other types of communication systems. 

The network part of the communication system includes a base transceiver station 
14 with which the mobile station 12 communicates by way of radio channels defined 
upon a radio air interface extending therebetween. Here, arrows 16 and 18 are 
representative of a forward link and a reverse link, respectively, upon which forward 

20 link and reverse link channels are defined. Data originated at the mobile station 12 for 
communication to the network part of the communication system is converted into a 
form to permit its communication upon reverse link channels defined upon the reverse 
link 18. And, analogously, data originated at the network part for communication to the 
mobile station is converted into a form to permit its communication upon forward link 

25 channels defined upon the forward link 16 to the mobile station. 

The network part of the communication system includes further functional 
entities, here including a base station controller (BSC) 22 and a mobile switching center 
(MSC) 24. The base station controller is coupled to groups of base transceiver stations, 
of which the base transceiver station is representative. The base station controller 

30 operates, amongst other things, to control operation of the base transceiver stations to 
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which the base station controller is coupled. The mobile switching center, in turn, is 
coupled to groups of base station controllers, of which the base station controller 22 is 
representative, and forms a switching center, or otherwise forms a gateway to external 
networks, here a PSTN (public switched telephonic network) and a packet data network 
5 (PDN) 28. A correspondent node (GN) 32 is coupled to both of the networks. The 
correspondent node is representative of a data source or a data sink that forms an 
ultimate source or destination of data that is communicated with the mobile station 12. 

The mobile station and the base transceiver station each include radio transceiver 
circuitry permitting both transmission and reception of data by way of the radio air 

10 interface formed therebetween. The base transceiver station, accordingly, includes a 
transmit part 36 and a receive part 38. And, the mobile station includes a transmit part 
42 and a receive part 44. Data to be communicated by the base transceiver station to the 
mobile station is communicated by the transmit part 36 and detected by the receive part 
44 of the mobile station. And, the data that is to be communicated by the mobile station 

15 is communicated by the transmit part 42 upon the radio air interface to be detected by 
the receive part 38. Operation of an embodiment of the present invention shall be 
described with respect to the generation and communication of data upon the forward 
link by the transmit part 36 of the base transceiver station for delivery to the receive part 
of the mobile station. Operation of the transmit part 42 and receive part of the mobile 

20 station and base transceiver station operate analogously. 

Data that is to be communicated by the network part to the mobile station is 
provided to the base transceiver station by way of the line 48. The data is applied to a 
source encoder 52. The source encodes the data and generates source encoded data on 
the line 54. Representations of the source encoded data is provided to a modem 56, 

25 here formed of a channel encoder 58 and a modulator/mapper 62. The data applied to 
the modem is channel encoded by the channel encoder and modulated by the modulator 
62. Once modulated, the data is provided to an antenna transducer 64. When the base 
transceiver station is of a multiple antenna formation, the modulator also forms mapping 
functions to map the data to selected ones of the antennas 64. 
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Once transduced into electromagnetic form, the data is communicated upon one 
or more forward link channels to be detected by the receive part 44 of the mobile 
station. The receive part of the mobile station includes functional elements 
corresponding to, but operating generally reverse to that of, the transmit part of the base 
5 transceiver station. More particularly, the receive part of the mobile station includes a 
modem 68 that includes a demodulator 72 and a channel decoder 74. The demodulator 
demodulates data detected at the antenna 76 of the mobile station. The channel decoder 
channel decodes the demodulated data applied thereto. A source decoder 78 also forms 
part of the receive part of the mobile station. 

10 In the exemplary implementation, digital communication techniques are utilized, 

and a communication service effectuated pursuant to operation of the communication 
system is effectuated through the communication of sequences of data symbols that 
together form the data that is to be communicated between communication stations of 
the communication system. Mathematical operations must be performed upon 

15 successive symbols of the data sequence. For instance, at the modem devices 56 and 68, 
determinations of Euclidean norms associated with the successive data symbols must be 
made. As noted previously, Euclidean norm calculations are computationally intensive, 
sometimes prohibitively so. Manners by which to approximate, or otherwise estimate, 
the Euclidean norm, are therefore needed. Existing schemes are, however, practically 

20 unavailable when higher dimensional Euclidean norms are to be approximated. And, in 
two dimensions, improved accuracy of the approximation would be beneficial. 

Accordingly, apparatus, shown generally at 82, of an embodiment of the present 
invention is provided by which to approximate Euclidean norms in any of N dimensions. 
Two dimensional, three dimensional, etc. Euclidean norms are approximated through 

25 operation of the apparatus 82. At higher dimensions, approximations are performed 
iteratively, taking advantage of prior calculations to form successive two dimensional 
approximations, based upon prior calculated approximations during the iterative 
process. 

The apparatus 82 is formed of functional entities, here shown in functional form, 
30 but implementable in any desired manner. In one implementation, for instance, the 
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apparatus 82 forms an ASIC. In another implementation, the apparatus is implemented 
as algorithms executable by processing circuitry. Other implementations of the 
' apparatus are analogously possible. Here, the apparatus is formed of a transformer 84 
an estimator 86. The transformer 84 is provided, here on the line 88, with indications of 
5 component portions of a data symbol that is to be operated upon to obtain an 
approximation of the Euclidean norm thereof. The transformer transforms the 
component portions into transformed values, and transformed values are provided to the 
estimator 86. The estimator estimates a Euclidean norm based upon the transformed 
values provided thereto. A feedback loop 92 is also represented in the Figure. When a 

10 higher level dimension, i.e., greater than two, the feedback loop is utilized to provide 
indications associated with the estimated Euclidean norm in a subsequent iteration of 
the approximation process. When only a two dimensional Euclidean norm is to be 
formed, iterative operation of the apparatus is not required, and the loop 92 is not 
utilized. Analogous structure of the apparatus 82 is also embodied at the mobile station 

15 12, also operable to approximate N dimensional Euclidean norms, as necessary. 

When QPSK modulation techniques are utilized, two dimensional quadrature 
signal pairs are applied to the apparatus 82 and the approximations of the Euclidean 
norm thereof are made in a single pass. When higher level modulation techniques are 
utilized, iterative operation of the apparatus is performed. 

20 Mathematical analysis and operation of the apparatus follows. 

In the exemplary implementation, the design is to approximate the magnitude of a 

quadrature signal pair (i.e., approximate: R = ^I 2 +Q 2 , where / and Q represent the 

quadrature signal pair). First, since only the magnitude is needed, a the following 
simple transformation to rotate the complex point z = / + jQ such that its argument lies 
25 between 0 and tt/4 is made 

x = max(j/|,|0|) 

y = minj/|,$). (1) 

30 And therefore 
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R = Jx 2 +y 2 



0<,9<kI4 



(2) 



10 



15 



20 



25 



Since a fixed magnitude lies on an arch in the 2-D space the following approximation is 
made 



R = ax + by = R(a cos 6 + b sin 0) 
Thus the relative error is given by 



(3) 



R — R 

e{0) = = l-a cos 0 - b sin 0 

R 



Next the arch described in (3) is then divided into two regions 



(4) 



region I : 
region II : 



O<0<0 o 
0 o <0<tt/4 



(5) 



This division into two regions is used to improve the approximation by optimizing the 
solution based on the location of the complex point. Thus two solutions are formed and 
depending on the location the appropriate equation is chosen. 

Figure 2 illustrates the two-region approach showing both the actual magnitude 
and the approximation. Two regions, region I 96 and region II 98, are represented. 
Since the two arcs intersect at two points in each region there are four values where the 
error is zero. Additionally, it is clear that error achieves its peak amplitude at three 
points in each region, the mid and end points. The equiripple error criteria is meet in 
the two regions by setting 



region!: )e(0)\ = ^{0^ )| = \e(0 o )| 
region II : \e(0 o )| = \e{0^ n )| = |e(*/4)| 

where, as can be seen in Figure 2, for region I: \e{0 o \ = 



2) 



(6) 



= |e(;r/ 4 )| an d for region II: 



<e 0 } = 



=K*/4 



30 
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Using the expressions in (4) and (6) the following set of equations are formed 
Region I: 



Region II: 



e(0)| = l-a 
- \e(0 o /2)| = 1 - a cos(0 o /2) - b sinfo /2) 
le^^l-acos^-ftsin^ 

|e(^ 0 )| = l-acos^-6sin^ 



(7) 



icjA-9 0 



~\- a cos 



-b sin 



(8) 



|e(^/4)| = l - a cos(^/4) - b sin(>r/4) 



First region I is considered. Defining, for clarity, e, = |<O)|=je(0 o /2)| = |e(0 o ) reordering 

10 terms and using the double angle identities: cos2^ = 2cos 2 A-l and sin 2^ = 2 cos ^ sin ^ , 
equation (7) can be rewritten as 



15 



20 



Region I: 

bsin(0 o /2) + acos(0 o /2) - e 7 -l = 0 

26sin(0 o /2)cos(0 o /2)+a{2cos 2 (0 o /2)-l}+ e, -1 = 0 
0 + a + e, -1 = 0 

Next if line 1 of (9) is multiplied by -2cos(<9 0 /2) and then added to line 2 of (9) 

fcsin(0 o /2) + acos(0 Q /2) - e, - 1 =0 

0 -a + e 7 (2cos(0 o /2)+l)+(2cos(0 o /2)-l) = O 

0+ a + e, - 1 =0 



Continuing, line 2 is added to line 3 of (10) 



6sin(0 o /2) + acos(0 o /2) - e, - 1 =0 

0 -a +' e 7 (2cos(0 o /2)+l)+(2cos(0 o /2)-l)=O 

0 + 0 + - 2e 7 (cos(e 0 /2)+l) -2(cos((9 0 /2)-l) = 0 



25 Solving line 3 of (1 1) gives 
e,=tan 2 {0 o /4} 



(9) 



(10) 



(11) 



(12) 
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where again the double angle identity was used. 

Using (4) and (6), as was done above to determine j| , both a and b can also be 

found. Since the math is tedious the detailed solution for a and b will also be omitted. 
The solution, however, provides the following 



Region I: 



Region II: 



a = 



l + secfo/2) 
6 = 2tan(0 o /4) 

K0n*x,/)| = tan 2 (0 o /4) 



(13) 



a = - 



s/2(l-V2sinfl g ) 



2sini|^-e o ]+cos[^+^ 



b = 





j + cos 











(14) 



^ m ax,/)| = tan 2 i^ + ^ 0 j 



The maximum error percentage for Region I and Region II are given in (13) and 
(14) respectively. Figure 3 plots the two error curves 100 and 102 as a function of the 
angle of separation. Additionally, equating the two error expressions in (1 3) and (14) 
provides the optimal selection of the angle of separation: 
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|e(*max,/)|=K<W)| 
tan 2 (^/4) = tan 2 l^ + ^ 



0 



(15) 



As seen in Figure 3, the point at which both regions have equal maximum errors 
is 0.393 radians as was shown in (15). Although this angle indeed provides the best 
solution it does so with some implementation complexity. 

To determine in which of the two regions the argument lies requires a division, 
and the number of bits in that division must be reasonable large. Other techniques can 
be used to determine the region, which are more conducive to an ASIC 
implementation, but which do so at some cost in performance. One approach that 
readily allows this reduction is to make 0 o z function of tan _1 () . This assignment is 
clearly advantageous when considering that the desire is to approximate the magnitude 
of a quadrature signal. Two likely candidates are 6 0 = tan" 1 (^) and 0 O =tan _1 (|) these 

require only shift operations of the quadrature components in order to determine the 
region. 

This in itself does not completely address the entire complexity issue. Solving 
for the coefficients a and b using either the optimal, or one of the two aforementioned 
approximations, for 0 O requires a non-multiple-of-two multiplication. Therefore 
quantization of a and b is required. 

The following table shows the mathematical process, and its hardware 
implementations, by which to approximate a 2D-Euclidean norm. 

Fast Algorithm Details 



Fast Algorithm 



Fast Algorithm - Detailed 
Hardware implementation 
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x = max! 



>> = min(j/|,|2|). 
if(x > Ay) 

R = x 
else 

R =—x+—y 
8 2 



//(|/|>=|e|) 

x = \l\ 

. Hel 

else 

X = \I\ 

if\x>{y«2)} 

R = x 
else 

R = x-(x »3) + (y »i) 



Pursuant to another implementation, the following tables show the mathematical 
processes, and hardware implementations, by which to approximate two dimensional 
Euclidean norms. 
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Algorithm 


Algorithm - Detailed Hardware 




Implementation 


* = maxj/|,|e|) 
y = min(\l\,\$). . 

if(x>=2y) 

A 63 1 

R = — x+—y 
64 4' 

else 

n 7 1 

R=—x+—y 
8 l" 


//(|/|>=|e|) 

X = \I\ 

else 

*/{*>(>;«l)} 

# = x-(;c»6) + (j>»2) 

£ = jc-(*»3) + (j;»1) 




Algorithm 


Algorithm - Detailed Hardware 




Implementation 


* = max(j/|,|e|) 

j = min(j/|,|e|). 

z/(x>=2j) 

- 31 1 
R = — x+—y 
32 4 

S 7 1 
8 2 y 


//(|/|>=|e|) 

X = \I\ 

else 

H'l 

*/{x>0>«l)} 

# = *-(;<: » 5) + (y»2) 
else 

R = x-(x»3) + (y»\) 



In another implementation, particularly amenable for use when the apparatus is 
DSP-, or otherwise software-, implemented, the procedures set forth in the following 
table are carried out. 
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Algorithm - Detailed 
Hardware Implementation 



x = max(j/|,|e|) 
y = min(j/|,|0|). 

if{x >= 0.3927;;) 

2 / x 

R = , , v x + 2tan(fl 0 /4)v 

else 

V2"(l -TlsinflJ 



2sinl^-^j + cos^ + ^j 



■J2~(j2cos0 o -l) 



2sinl --0„ +cos 



else 
x = l 7 l 

>= 0.3927j>) 
£ = .9903* + .197>> 

= . 8395a: + . 561 Oy 



15 



Approximation of the Euclidean norm is also possible in three dimensions. In the 
exemplary implementation, the implementation set forth above is extended to three 
dimensions.. 

For clarity, the development of the N-dimensional solution will be presented 
using a 3-D example. This will then be generalized for higher dimensions. The problem 
of approximating the Euclidean norm is depicted in Figure 4 where the norm is 



represented by 



OP 



The point, P is located at [x,y,z] in rectangular coordinates. 



10 Mapping to spherical coordinates is done as: 

* = i?sin^cos0, y = /?sin^sin0, z = Rcos0 (16) 
Again, as in the 2-D problem, the task is to estimate the Euclidean norm 



O<0<;r/4, 0<^<^/2. 



(17) 



Where now three arguments are used for the three dimensions. If the same approach as 
in the 2-D problem an estimate may be formed as 



20 



R = ax + by + cz = R(a cos0sin^ + Z?sin0sin^ + ccos^) 



(18) 
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Therefore the relative error is given by 
R — R 

e(0) = —— = l-acos0sin0-bsm0sm0-ccos0 (19) 
5 The equiripple error criteria is meet in the two regions by setting 
region / : |e(0,0)| = |e(^, , )| = \e(0 o , ^ )| 

(20) 

region//: |^,#o)|=j^W.^)| = K*/W4)| 

The solution is attained by solving (19) and (20) for the unknowns a,b and c with 
10 the assignment of 0 O =<fi 0 = ;r/8as was found previously for the 2-D problem. This 

solution is relatively complex and it is questionable whether the approximation approach 
is less complex then the actual solution. 

The approach to solve the 3-D problem is vastly more complicated then that of 
the 2-D problem. Thus a novel approach was used whereby the 3-D problem was 
15 reduced to two 2-D problems and then solved iteratively. The solution is as follows: 
First a set of 2-D components are used to find the magnitude in 2-D space using one of 
the approaches outlined previously. Next, the resultant magnitude and remaining 
unused 3-D component are used to form a second 2-D set. Again using one of the 
approaches outline in Section 2 the magnitude is computed to provide a final 
20 approximation to the original 3-D Euclidean norm. 

Stated another way, first a 2-D magnitude is found using the points (x,y) 9 then the 
x-axis (arbitrary selection) is rotated by0 to the resultant vectors position (o y P') creating 
the x r axis. Next a second 2-D problem is solved using the x' and z axis. 

The following steps outline the algorithm used to find the Euclidean norm for a 
25 3-dimension problem given (X f Y y Z); 

1.) Find: 

x = max 



H\y\) 

30 y = min{4|y|). 
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2. ) Approximate the 2-D magnitude, R , using one of the approaches outlined 
previously. 

3. ) Next, rotate the x-axis to form thex' axis and define: 

5 

x = max|ff|,|z|) 
: >/=mm(tf|,|Z|). ' 

4. ) Next approximate this second 2-D magnitude using again one of the 

10 approaches from Section 2 for R attaining the final 3-D Euclidean norm approximation. 
Figure 5 illustrates a method, shown generally at 132, representative of the 
method of operation of an embodiment of the present invention. The method forms an 
approximation of the Euclidean norm of a data symbol defined in terms of a first 
component portion and at least a second component portion. 

15 First, and as indicated by the block 134, indications of the first and at least 

second component portions of the data symbol are detected. Then, and as indicated by 
the block 136, the first component portion is transformed into a first transformed value 
and the second component portion is transformed into a second transformed value. The 
first and the second transformed values are selected such that the geometric argument 

20 defined therefrom is within a selected angular range. Then, and as indicated by the 
block 138, the Euclidean norm of the data symbol is estimated. The Euclidean norm is 
estimated to be the summation of the first product value formed from the first 
transformed value and a second product value formed from the second transformed 
value. 

25 Higher level Euclidean norm approximations are performed by an iterative 

procedure in which two dimensional Euclidean norm approximations are successively 
made to form the approximation of the N dimensional Euclidean norm. Accurate 
approximations are obtained without necessitating excessive computations. Thereby, 
approximations of Euclidean norms at higher dimensions are possible at high rates, such 

30 as in modem devices operable to form parts of communication stations of a cellular 
communication system that provides high speed data services. 
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The previous descriptions are of preferred examples for implementing the 
invention, and the scope of the invention should not necessarily be limited by this 
description. The scope of the present invention is defined by the following claims. 
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